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Abstract 

Using exact enumerations of self-avoiding walks (SAWs) we compute the inhomogeneous pressure 
exerted by a two-dimensional end-grafted polymer on the grafting line which limits a semi-infinite 
square lattice. The results for SAWs show that the asymptotic decay of the pressure as a function of 
the distance to the grafting point follows a power-law with an exponent similar to that of gaussian 
chains and is, in this sense, independent of excluded volume effects. 

PACS numbers: 05.50.+q,36.20.Ey 
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I. INTRODUCTION 



Imaging and manipulating matter at sub-micron length scales has been the cornerstone 
of nano-sciences development [lj. In Soft Matter systems, including those of biological 
relevance, the cohesive energies being only barely larger than the thermal energy ksT, forces 
as small as a pico- Newton exerted over a nanometer length scale might be significant enough 
to induce structural changes. Examples can be found in the stretching of DNA molecules by 
optical traps [2j, on the behavior of colloidal solutions under external fields [3J and on the 
deformations of self-assembled bilayers |4j] to name just a few. Thus, in Soft Matter, when 
one exerts a localized force over a small area, precise control of the acting force requires not 
only a prescribed value of the total applied force but, more importantly, a precise pressure 
distribution in the contact area. 

The microscopic nature of pressure has been understood since the seminal work of 
Bernoulli two and a half centuries ago: in a container, momentum is transferred by collisions 
from the moving particles to the walls 5] . When the particle concentration is homogeneous 
so is the pressure. Strategies for localizing the pressure over a nanometer area thus requires 
the generation of strong concentration inhomogeneities, at equivalently small scales. Bickel 



et al. la iTfl and Breidnich et al. [8| have recently realized that such inhomogeneities are in- 
trinsic to entropic systems of connected particles such as polymer chains, and have computed 
the inhomogeneous pressure associated with end-grafted polymer chains within available an- 
alytical theories for ideal chains. Their results show that the polymer produces a local field 
of pressure on the grafting surface, with the interaction being strong at the anchoring point 
and vanishing far enough from it. Scaling arguments were also put forward in {t| to discuss 
the more relevant case of real polymer chains, where excluded volume interactions between 
the different monomers need to be taken into account. These arguments suggest that the 
functional variation of pressure with distance from the grafting point should be the same in 
chains with or without excluded volume interactions, albeit with different prefactors. 

In this paper we compute the inhomogenous pressure applied to a wall by an end-grafted 
polymer with excluded volume interactions, modeled by selfavoiding walks (SAWs) on the 
square lattice. In Fig. [TJ we illustrate our model with a wall located at x — 0. The wall is 
neutral, in the sense that the statistical weight of a monomer placed on the wall is equal to 
the weight of a monomer in the bulk. The length of a step of the walk is equal to the lattice 



constant a, and we use this as the length unit. The model is athermal, that is, all allowed 
configurations of a SAW have the same energy. 
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FIG. 1. A SAW grafted at the origin x = y = to a wall placed on the y axis. If the vertex on the 
wall at (0, 1) is not excluded, the only possibility for the next step would be towards this vertex. 
If this vertex is excluded, the SAW will end at the final point (1,1). 



The canonical partition function of walks with n steps (Z n ) is equal to the number of 
SAWs starting at the origin and restricted to the half-plane x > 0, called ch^ in [9|. The 
Helmholtz free energy is given by F n = -k B T\nc { n ] . We can estimate the pressure exerted 
by the SAW at a point (0,r) on the wall by excluding this vertex from the lattice. The 
excluded vertex is represented as a hatched square in Fig. [TJat r = 1. The pressure P n (r) 
exerted at this point is then related to the change in the free energy when the vertex is 
excluded, P n a 2 = —AF n . If we call Cn\r) the number of n step SAWs with the vertex at 
(0,r) excluded, the dimensionless reduced pressure may be written as 

P ra (r)a 2 d l \r) 

KB-L Cn 

Of course we are interested in the thermodynamic limit p(r) = lim„_ s>00 p n (r), so the enu- 
meration data must be extrapolated to the infinite length limit. It is worth noting that the 
density of monomers at the vertex (0,r) is given by p(r) = 1 — lim^oo Cn ^ (r)/cn \ so that 



p(r) = — ln[l — p(r)}. 
3 



(2) 



The exact enumerations allow us to obtain precise estimates of the pressure exerted by 
SAW's at small distances of the grafting point, and we find, rather surprisingly, that the 
asymptotic form of this pressure is well reproduced even for these small values of r. In 
section |TT] we give some details of the computational enumeration procedure. In section 
IIHI the enumeration data are analyzed and estimates for the pressure as a function of the 
distance to the grafting point are presented. Final discussions and conclusions may be found 
in section IIVI 



II. EXACT ENUMERATIONS 



The algorithm we use to enumerate SAWs on the square lattice builds on the pioneering 

n 

work of Enting [10|] who enumerated square lattice self-avoiding polygons using the finite 
lattice method. More specifically our algorithm is based in large part on the one devised 



by Conway, Enting and Guttmann 



ll| for the enumeration of SAWs. The details of our 
algorithm can be found in 12| . Below we shall only briefly outline the basics of the algorithm 
and describe the changes made for the particular problem studied in this work. 

The first terms in the series for the SAWs generating function can be calculated using 
transfer matrix techniques to count the number of SAWs in rectangles W vertices wide and L 
vertices long. Any SAW spanning such a rectangle has length at least W + L — 2. By adding 
the contributions from all rectangles of width W < N + 1 and length W < L < N — W + 1 
the number of SAW is obtained correctly up to length N. 

The generating function for rectangles with fixed width W are calculated using transfer 
matrix (TM) techniques. The most efficient implementation of the TM algorithm generally 
involves bisecting the finite lattice with a boundary (this is just a line in the case of rect- 
angles) and moving the boundary in such a way as to build up the lattice vertex by vertex 
as illustrated in Fig. [2J If we draw a SAW and then cut it by a line we observe that the 
partial SAW to the left of this line consists of a number of loops connecting two edges (we 
shall refer to these as loop ends) in the intersection, and pieces which are connected to only 
one edge (we call these free ends). The other end of the free piece is either the start-point 
or the end-point of the SAW so there are at most two free ends. 

Each end of a loop is assigned one of two labels depending on whether it is the lower 
end or the upper end of a loop. Each configuration along the boundary line can thus be 
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FIG. 2. A snapshot of the boundary line (dashed line) during the transfer matrix (TM) calculation 
on a strip of width 7 with r = 3. The filled circle indicates the grafted start-point of the SAW and 
the shaded box the excluded vertex. SAWs are enumerated by successive moves of the kink in the 
boundary line, as exemplified by the position given by the dotted line, so that one vertex and two 
edges at a time are added to the strip. To the left of the boundary line we have drawn an example 
of a partially completed SAW. 

represented by a set of edge states {crj}, where 

empty edge, 

1 lower loop-end, 

<T; < (3) 

2 upper loop-end. 

3 free end. 

If we read from the bottom to the top, the configuration or signature S along the intersection 
of the partial SAW in Fig. [2] is S = {031212120}. Since crossings aren't permitted this 
encoding uniquely describes which loop ends are connected. 

The sum over all contributing graphs is calculated as the boundary is moved through the 
lattice. For each configuration of occupied or empty edges along the intersection we maintain 
a generating function G$ for partial walks with signature S. In exact enumeration studies 
such as this G$ is a truncated polynomial G${x) where x is conjugate to the number of 
steps. In a TM update each source signature S (before the boundary is moved) gives rise to 
a few new target signatures S' (after the move of the boundary line) and m = 0, 1 or 2 new 
edges are inserted leading to the update Gs>(x) = Gs>(x) + x m Gs(x). Once a signature S 



has been processed it can be discarded. The calculations were done using integer arithmetic 
modulo several prime numbers with the full integer coefficients reconstructed at the end 
using the Chinese remainder theorem. 

Some changes to the algorithm described in 12] are required in order to enumerate the 



restricted SAW we study here. Grafting the SAW to the wall can be achieved by forcing the 
SAW to have a free end (the start-point) on the top side of the rectangle. In enumerations 
of unrestricted SAW one can use symmetry to restrict the TM calculations to rectangles 
with W < N/2 + 1 and L > W by counting contributions for rectangles with L > W twice. 
The grafting of the start-point to the wall breaks the symmetry and we have to consider all 
rectangles with W < N + 1. Clearly the number of configurations one must consider grows 
with W. Hence one wants to minimize the length of the boundary line. To achieve this the 
TM calculation on the set of rectangles is broken into two sub-sets with L >W and L < W, 
respectively. The calculations for the sub-set with L > W is done as outlined above. In 
the calculations for the sub-set with L < W the boundary line is chosen to be horizontal 
(rather than vertical) so it cuts across at most L + 1 edges. Alternatively, one may view the 
calculation for the second sub-set as a TM algorithm for SAW with its start-point on the 
left-most border of the rectangle. 

Exclusion of the vertex at distance r from the starting point of the SAW is achieved by 
blocking this vertex so the walk can't visit the vertex. The actual calculation can be done 
in at least two ways. One can simply specify the position of the starting point (and r) on 
the upper/left border and sum over all possible positions. This means doing calculations 
for a given width W many times; once for each position of the starting point of the SAW. 
Alternatively one can introduce 'memory' into the TM algorithm. Specifically once we have 
created a configuration which inserts the first free end we 'remember' that it did so. We can 
flag that the free end has been inserted by adding a ghost edge to the configuration initially 
in state 0. Once the first free end is inserted the state of the ghost edge is changed to 1. 
In the next sweep the state of the ghost edge is incremented by 1. When the state of the 
ghost edge has reached the value r the vertex on the top border is blocked. The problem 
with the first approach is that we need to do many calculations for any given rectangle. 
The problem with the second approach is that we need to keep r + 1 copies of most TM 
configurations thus using substantially more memory. The choice will be a matter of whether 
the major computational bottle-neck is CPU time or memory. For this study we used the 



first approach. 

In more detail the TM algorithm for the case L >W works as follows. A SAW has two 
free ends and in the TM algorithm the first free end is forced to be at the top at a distance 
k from from the left border (this is the starting point of the SAW). We then add a further 
r — 1 columns; in the next column the top vertex is forced to be empty. After this further 
columns are added up to a maximum length of L m = N — W + 1. This calculation is then 
repeated for k — to L m thus enumerating all possible SAWs spanning rectangles of width 
exactly W and length L > W. A similar calculation is then done with the SAW grafted to 
the left border and in each case repeated for all W < N/2. 

The calculation above enumerates almost all possible SAWs. However, we have missed 
those SAWs with two free ends in the top border where the end-point precedes the starting- 
point. That is there is a free end in the top border at a distance > r prior to the excluded 
vertex. We need to count such SAWs separately. The required changes to the algorithm are 
quite straight-forward and will not be detailed here. 

We calculated the number of SAWs up to length n = 59 for the unrestricted case and for 
an excluded vertex with r = 1, 2, 3, 4, 5, 10, 20. In each case the calculation was performed 
in parallel using up to 8 processors, a maximum of some 16GB of memory and using a total 
of under 2000 CPU hours (see 12| for details of the parallel algorithm). We needed 3 primes 
to represent each series correctly and the calculations for all the primes were done in a single 
run. 



III. ANALYSIS AND RESULTS 

In tables UJ and [Til we have listed the results for the enumerations of self-avoiding walks 
without additional restrictions, ch \ and walks which are not allowed to occupy the vertex 
(0, 1) of the wall, c£\l). If we calculate the pressures directly, we notice a parity effect, as 
seen in the results presented in Fig. [3] This effect is related to an unphysical singularity 
in the generating function of the counts c£ , G(x) = 2~2^=o c ^ x n . Besides the physical 
singularity at x = x c = l//i, where [i is the connective constant, there is another singularity 



at x = — 1//J, This point will be discussed in more detail below, and more precise 
estimates for the pressures at several distances from the grafting point will be provided. 




FIG. 3. Pressure p n ( r ) for r = 1, calculated with the enumeration data for ch and c~ n (1) using 
expression ([!]). 



A. Critical points and exponents 

The critical behaviour of a polymer grafted to a surface is well established [13[. It has 
)een proved that the connective constant of grafted walks equals that of unrestricted walks 
Ijj]. The associated generating function has a dominant singularity at x = x c = 



G{x)=Y, C n )xn ~ A ^-^y 



71 



(4) 



where 71 = 61/64 is a known Q| critical exponent. Besides the physical singularity 



there is another singularity at x 



x. 



We have analysed the series using differential approximants 18J. We calculate many 
individual approximants and obtain estimates for the critical points and exponents by an 



averaging procedure described in chapter 8 of reference 



19j. Here and elsewhere uncertainties 



on estimates from differential apprimants was obtained from the spread among the various 



191 ] . The results for unrestricted grafted SAWs are listed in 

1, 2, 5 and 10.. From 



approximants as detailed in 

Table II I II under r = 0. We also list estimates for the cases r 



these estimates it is clear that all the series have the same critical behaviour. That is a 
dominant singularity at x = x c with exponent —71 = —61/64 and a non-physical singularity 
at with a critical exponent consistent with the exact value 7_ = 3/2. 

The critical behavior can be established more rigorously from a simple combinatorial 
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TABLE I. Number of walks in the half-plane o, 



n 




n 


J 1 ) 


n 


(X) 


1 


3 


21 


681552747 


41 


176707555110156095 


2 


7 


22 


1793492411 


42 


465629874801142259 


3 


19 


23 


4725856129 


43 


1227318029107006037 


4 


49 


24 


12439233695 


44 


3234212894649555857 


5 


131 


25 


32778031159 


45 


8525055738741918835 


6 


339 


26 


86295460555 


46 


22466322857670716727 


7 


899 


27 


227399388019 


47 


59220537922987286933 


8 


2345 


28 


598784536563 


48 


156073168859898607113 


9 


6199 


29 


1577923781445 


49 


411414632591966686887 


10 


16225 


30 


4155578176581 


50 


1084313600069268939547 


11 


42811 


31 


10951205039221 


51 


2858360190045390998925 


12 


112285 


32 


28844438356929 


52 


7533725151809823220637 


13 


296051 


33 


76016486583763 


53 


19860118923927104821817 


14 


777411 


34 


200242023748929 


54 


52346889766180530489735 


15 


2049025 


35 


527735162655901 


55 


137997896899080793506959 


16 


5384855 


36 


1390287671021273 


56 


363744527134008049572583 


17 


14190509 


37 


3664208598233159 


57 


958930393586321187515995 


18 


37313977 


38 


9653950752700371 


58 


2527696511232818406275131 


19 


98324565 


39 


25444550692827111 


59 


6663833305674862002802763 


20 


258654441 


40 


67042749110884297 







argument. The number of walks Cn\r) with the point at (0,r) excluded is clearly less than 
the number of unrestricted walks Cn\ On the other hand if we attach a single vertical step 
to the grafting point of an unrestricted walk we get a walk which does not touch the surface 
at all and hence these walks are a subset of Cn\r). This establishes the inequality 

c£>i < cf>(r) < c« (5) 

and hence shows that up to amplitudes the asymptotic behaviors of these sequences are 
identical. 
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TABLE II. Number of restricted walks in the half-plane c n (1). 



n 


<V(1) 


n 


(l),^ 

Cn {1} 


n 


da' {I) 


1 


2 


21 


A i~\ A O C~\ f~\ O 

484553893 


41 


125845983216200025 


2 


5 


22 


1277403184 


42 


331741159147128245 


3 


13 


23 


3361118347 


43 


874112388226242422 


4 


35 


24 


8860136085 


44 


2304278197456842952 


5 


91 


25 


23319106552 


45 


6071977423574762560 


6 


242 


26 


61468398004 


46 


16006835327039914244 


7 


630 


27 


161814936995 


47 


42181825940070651834 


8 


1672 


28 


426530787110 


48 


111200914189945767681 


9 


4369 


29 


1123043680259 


49 


293056004233059019257 


10 


11558 


30 


2960232320818 


50 


772575890795109134325 


11 


30275 


31 


7795418415398 


51 


2036121996024316003415 


12 


79967 


32 


20548006324647 


52 


5367866589569286706072 


13 


209779 


33 


54117914172220 


53 


14147607361624429924807 


14 


v t* o r* o A 

553634 


34 


142651034798697 


54 


37298221266819312654286 


15 


1453801 


35 


375747632401071 


55 


98307470253293931954939 


1 p. 

ID 


ooo4o / O 


'ip. 

oO 


QQD/1 ^fi^n7ni 1 0.9Q 


KP. 
00 


zoy i / ooUoozuzoi izzy < 4Zou 


17 


10077384 


37 


2609158017850105 


57 


683144867659867533730505 


18 


26574366 


38 


6877742334133961 


58 


1801074652042354959971779 


19 


69870615 


39 


18119629209950641 


59 


4747450605648675761162683 


20 


184216886 


40 


47764129557587369 







B. Pressure 

Having established the critical behaviour of the series we can now turn to the determi- 
nation of the pressure exerted by the polymer on the surface. Since all the series have the 
same dominant critical behaviour it follows from (CQ) that the pressure is given by the ratio 
of the critical amplitudes. 

One way of estimating the amplitudes is by a direct fit to an assumed asymptotic form. 
Here we assume that the asymptotic behaviour of our series is similar to that of unrestricted 
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TABLE III. Estimates of the critical points and exponents for SAWs with an excluded vertex a 
distance r from the origin (r = is the unrestricted case). The estimates were obtained from third 
order approximants with L being the degree of the inhomogenous polynomial. 



r L x c 

0.379052260(64) 

4 0.379052241(20) 

8 0.379052243(14) 



7 

0.953097(70) 
0.953072(17) 
0.953071(15) 



-0.3790526(38) 
-0.3790492(30) 
-0.3790498(21) 



7- 

1.5002(19) 
1.5023(13) 
1.5016(12) 



1 0.3790522582(30) 
1 4 0.3790522575(38) 
1 8 0.379052257(11) 



0.9530884(24) 
0.9530879(30) 
0.953090(14) 



-0.3790425(97) 

-0.379030(26) 

-0.379058(16) 



1.5074(74) 

1.523(29) 

1.4988(69) 



2 0.379052292(16) 
2 4 0.379052276(12) 
2 8 0.379052306(26) 



0.953123(13) 

0.9531115(97) 

0.953135(20) 



-0.3790511(33) 
-0.3790478(89) 
-0.379057(21) 



1.5011(24) 
1.5036(60) 
1.498(20) 



5 0.37905218(21) 
5 4 0.37905225(31) 
5 8 0.37905226(29) 



0.95304(17) 
0.95313(24) 
0.95313(25) 



-0.379114(61) 
-0.379099(40) 
-0.379074(31) 



1.457(37) 
1.467(29) 
1.482(20) 



10 0.3790483(12) 
10 4 0.3790493(40) 
10 8 0.3790508(22) 



0.9494(12) 
0.9503(32) 
0.9514(14) 



-0.379230(55) 
-0.379237(29) 
-0.379246(91) 



1.369(32) 
1.369(14) 
1.365(54) 



SAW 20] . The asymptotic analysis of 20j was very thorough and clearly established that 
the leading non-analytic correction-to-scaling exponent has the value 3/2 (there are also 
analytic, i.e., integer valued corrections to scaling). We repeated some of the steps in this 
analysis with the same result for the leading non- analytic correction-to-scaling exponent. 
Naturally there may be further non-analytic correction-to-scaling exponents with values 
> 3/2, but these would be impossible to detect numerically with any degree of certainty. 
So here we assume that the physical singularity has a leading correction-to-scaling exponent 
of 1 followed by further half-integer corrections while we assume only integer corrections at 
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FIG. 4. Estimates for the leading amplitudes obtained by fitting to the asymptotic form ([6]) plotted 
against 1/n while truncating the asymptotic expansion after 4 to 7 terms. 



the non-physical singularity. We thus fit the coefficients to the asymptotic form 



n 



71' 



1 (A(r) + J2 a A r )/ nj/2 ) +(-l) n n^- 1 J2 b k(r)/n k 

\ j=2 J k=0 



(6) 



In the fits we use the extremely accurate estimate /i = 2.63815853035(2) obtained from an 
analysis of the series for self-avoiding polygons on the square lattice 2lJ and the conjec- 
tured exact values 71 = 61/64 and 7_ = 3/2. That is we take a sub-sequence of terms 
{cW ( r )> c i-i( r )> • • • > c i-2m-i( r )}> P m § into tne formula above taking m terms from both the 
cij and bk sums, and solve the 2m linear equations to obtain estimates for the amplitudes. 
It is then advantageous to plot estimates for the leading amplitude A{r) against 1/n 
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for several values of m as done in Fig. HJ The behaviour of the estimates for the leading 
amplitudes shown in this figure supports that (jSJ) is a very good approximation to the true 
asymptotic form. In particular note that the slope becomes very flat as n is increased and 
decreases as the number of terms m included in the fit is increased. From these plots we 
estimate A = 1.124705(5), A(l) = 0.801625(5), A(2) = 0.97564(2) and A(5) = 1.09325(10), 
where the uncertainty is a conservative value chosen to include most of the extrapolations 
from Fig. HI 

The amplitude ratios A(r)/A, and hence the pressure, can also be estimated by direct 
extrapolation of the relevant quotient sequence, using a method due to Owczarek et al. 



221 ]: Given a sequence {a n } defined for n > 1, assumed to converge to a limit a ro with 
corrections of the form a n ~ a 0O (l + 6/n+...), we first construct a new sequence {p n } 
defined by p n = Ylm=i a m- We then analyse the corresponding generating function 

Estimates for o m and the parameter b can then be obtained from differential approximants, 
that is cioo is just the reciprocal of the first singularity on the positive real axis of P(x). 
In our case we study the sequence of ratios a n (r) = Cn\r)/ch\ which has the required 
asymptotic form. Using the same type of differential approximant method outlined above 
we find that A/A(l) = 1.4030218(5), which is entirely consistent with the estimate A/A(l) = 
1.403030(15) obtained using the amplitude estimates from the direct fitting procedure. 
Next we compare these results for the pressure with the ones for gaussian chains as 



expressed in equation (4) in jg]. That expression is for polymers in a three-dimensional 
half-space confined by a two-dimensional wall, and corresponds to finite values of the radius 
of gyration. If the expression is generalized to the <i-dimensional case and restricted to the 
limit of infinite chains, where the radius of gyration diverges, the result is: 

_ P G (r)a d _ T(d/2) 1 
PG{) k B T vr rf / 2 (r 2 + l) rf / 2 ' {) 

where we recall that r is dimensionless, measured in units of the lattice constant a. In table 
HVlwe have listed the estimated pressures for SAWs and the pressures obtained for Gaussian 
chains in d = 2, on the semi-infinite square lattice. In Fig. Owe have plotted the pressure for 
polymers modelled as SAWs and as Gaussian chains. In this figure the dashed line represents 
a decay in pressure with the same asymptotic form, ex l/(r 2 + 1), as the Gaussian chain but 
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r p(r) — SAWs p(r)-gaussian 

1 0.33863 0.15915 

2 0.14218 0.06366 

3 0.07334 0.03183 

4 0.04347 0.01872 

5 0.02844 0.01224 
10 0.00735 0.00315 

TABLE IV. Pressure at a distance r from the grafting point for SAWs and Gaussian chains. 

normalised so the curve passes through the SAWs data point for r = 10. Quite clearly the 
SAWs data is well represented by this form even for small distances r > 2. For r = 20 the 
SAWs data was indistinguishable from zero pressure. 



IV. FINAL DISCUSSIONS AND CONCLUSION 



Since our model is athermal and discrete, it is not really possible to compare our results 
with those obtained for the gaussian chain. However, as was already mentioned by Bickel et 



al. 



the excluded volume interactions should not change the scaling form of the pressure. 



Fig. EJ^b) clearly shows a 1/r 2 decay of the pressure, even for small distances. According 
to Bickel et al. [7], this similarity is due to the fact that the pressure and the monomer 
concentration in the vicinity of the wall are linearly related. On the other hand, it seems 
that the concentration is not affected by the molecular details or by the differences between 
chain models. In our case, despite the fact that p(r) and p(r) are related by a logarithmic 
relation, as shown in expression (T5]), we have for r > 1 a small concentration leading to a 
linear relation between those quantities. Actually, even for r ~ 2, we can observe a linear 
dependence, as shown in Fig. [6j 

Since the grafted chain is in mechanical equilibrium, the force J 7 applied to the walk at 
the grafting point, which is in the negative x direction in Fig. (TJ should be equal to the 
sum of the forces applied by the wall at other contact points, which are in the positive x 
direction. Thus, the dimensionless force is given by: 

° r=l 
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FIG. 5. (a) The pressure p(r) exerted by a polymer on a surface at a distance r from the grafting 
point. Data are for polymers modelled as SAWs or Gaussian chains. The dashed line are a 1/r 2 
fit. (b) Both data have the same 1/r 2 scaling form, even for values close to r = 2. The dashed 
lines are guide lines with slope equal to —2. 



For gaussian chains, integrating equation (J7J), we find fa = 1. For SAWs, we may estimate 
the force summing the results for r = 1,2, ... ,5 and obtaining the remaining contributions 
(r = 6, 7, . . . , oo) using the asymptotic result p(r) w A p /(r 2 + 1) where A p « 0.74235 was 
estimated using the result of p(r) for r = 10. The result of this calculation is fsAw ~ 
1.533, larger than the one for gaussian chains. As mentioned above, it does not seem 
straightforward to compare the two models, since a gaussian chain is a mass-spring model 
and therefore it is, unlike SAWs, not athermal. We may also mention that if p(r) for SAWs 
is extended to real values of r using a numerical interpolation procedure and the data for 
gaussian chains are rescaled so the areas below both curves are the same, the difference 
between the curves is quite small, the maximum being close to the origin and of order 
10 -3 . Due to the limited precision of the estimates for SAWs and to the expected small 

15 



0.4 



0,3 
On 0,2 
0,1 



"0 0,05 0,1 0,15 0,2 0,25 

P 

FIG. 6. Relation between pressure and concentration of monomers near to the wall at a distance 
r from the grafting point. For r > 1, a linear relationship is observed. 

dependency of the results on the interpolation procedure we will not present these results 
here, but we found that in general the rescaled results for the pressure of gaussian chains 
are larger than the pressures for SAWs at small values of r, but the inverse situation is 
found for larger distances. This net effect may be understood if we recall that the pressure 
is a monotonically growing function of the local density at the wall (Eq. and that the 
effect of the excluded volume interactions should be a slower decay of this density with the 
distance from the grafting point, as compared to approximations where this interaction is 
neglected. 

It is of some interest to obtain the total force applied to the chain at the grafting point for 
ideal chains, modeled by random walks on the semi-infinite square lattice. This force may 
be calculated considering the shift of the grafting point by one lattice unit in the positive x 
direction in Fig. HJ The change in free energy under this operation will be proportional to 
the force. This calculations should lead to the same result of the ones above, where the force 
was obtained summing over the pressures at all other sites of the wall besides the origin, 
since the total force applied on the chain has to vanish. 

Let us start by briefly reviewing the calculation of the number of random walks on a half- 
plane of the square lattice. If we call c n (p) the number of n-steps random walks on a square 
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lattice starting at the origin and ending at the point p — xi + yj, the number of RWs on 
the half-plane x > may be calculated by placing an absorbing wall at x — — 1, so that any 
walk reaching the wall is annihilated. This may be accomplished by using an image walker, 
starting at the reflection point of the origin with respect to the wall and ending at p. We 
will place the starting point of the random walk at (s, 0), where s = corresponds to walks 
starting at the origin. In this case the image walker starting point will be at po = — {s + 2)i, 
with distances measured in units of the lattice constant a. The number of walks confined to 



the x > half plane is given by 23] 



c«(p, s) = Cn{p) - c n (p+ (2 + s)i). (9) 

Since we are interested in the large n limit, we may use the gaussian approximation for the 
number of walks 

c n {p) = — expf-i^V (10) 



(11) 



nn \ n 

For the half-plane we get 



An 

c W(z s ) = — 



p\ 2 \ / |p+(2 + s)i| 2 



exp — exp 



n / \ n 



To obtain the total number of walks, we integrate this expression over the final point p 

c«( S )= / dx / dycW(p,s). (12) 

JO J-oo 

The result is 

4 n |.(2+s)/v^ 2 



/7T 



i 1} (*) 

for n ^> s, we have the asymptotic behavior 



e~ x dx, (13) 



$>{*) = V?^, (14) 



which has the expected scaling form (jlj), with exponent 7 = 1/2 and amplitude A = 2(s + 
1 ) / ^/tt. The change in free energy between the cases with s = and s — 1 is therefore given 
by — A: j bT' In 2, so that the force applied to the polymer by the wall at the grafting point will 
be fnw = In 2 m 0.6931, which is lower than the forces obtained for gaussian chains and 
estimated for SAWs. 

It should be mentioned that for SAWs the sum of the pressures corresponding to two 
distances p(rj) +p(rj) is always smaller (for finite — than — AF(ri, rj)/(/c B T), where 
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AF(ri,rj) is the change in free energy when both cells, at r, and r 3 - are excluded. In other 
words, an effective attractive interaction exists between the two excluded cells, so that the 
free energy decreases as the cells approach each other. This effect is due to walks in the 
unrestricted case which visit both excluded cells, and are therefore not counted in either 
Cn or Cn\r 2 ). The total force f' S Aw> resulting from the simultaneous exclusion of all 
cells besides the one at the grafting point r = 0, must thus be smaller than the force Jsaw 
defined in equation (JSJ). It is easy to find, since the number of SAWs with n steps So 1 in 
this case is given by dn ^ = 1 + G^lx, that for a given value of n the force at the grafting 
point will be f' nS Aw = ~ hi(c?n /cn )■ For large n, we get f' SAW — hx/i ~ 0.9701, smaller 
than fsAw — 1-533, as expected. 

Finally, we should also stress that although the pressure applied by the SAWs and by 
the gaussian chains display a similar power-law behavior, other possible walks on the lattice 
might lead to different results. Recently the pressure exerted by directed walks starting at 
the origin on the limiting line of a semi-infinite square lattice was obtained 24| . In the limit 
of large directed walks the asymptotic decay of the pressure with the distance to the grafting 
point also follows a power law, albeit with an exponent smaller than the one obtained here 
for SAWs and gaussian chains. 
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